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The effective coupling of QCD is measured by using the gauge configurations produced by the 
MILC collaboration in which the Kogut Susskind (KS) fermion is incorporated and by using that 
produced by the QCDOC collaboration in which the domain wall fermion (DWF) is incorporated. 
We fix the gauge to the Landau gauge and to the Coulomb gauge. The infrared effective coupling 
in the Coulomb gauge agrees with the recent extraction at JLab, but that in the Landau gauge 
shows infrared suppression. The suppression is expected to be due to the color anti-symmetric 
ghost propagator which in the unquenched configurations has stronger infrared singularity than 
the color diagonal ghost propagator. The Coulomb form factor in the infrared depends on the kind 
of the fermion incorporated in the system and the temperature. 

The quark has the effect of quenching randomness and the fluctuation of the color anti-symmetric 
ghost propagator is reduced in the unquenched configuration, and the Kugo-Ojima parameter c is 
closer to 1 in the unquenched configuration than in the quenched configuration. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 

t 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-Share Alike Licence. 



http://pos.sissa.it/ 



Roles of the quark field in the infrared lattice Coulomb gauge and Landau gauge QCD 



Sadataka Furui 



1. Introduction 

We fixed gauge configurations of the MILC collaboration!]]], §] and the QCDOC collaboration^] 
to Landau gauge and to Coulomb gauge and studied the infrared features of the lattice QCD by com- 
paring the gluon propagator, the ghost propagator, the effective coupling etc. In the case of finite 
temparature MILC f t , we assign the /3 = 5.65, 5.725 and 5.85 data correspond to T = 143, 172.5 and 
185 MeVp, and in the unit of critical temperature T c by T/T c = 1.02, 1.23 and 1.32, respectively. 
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2. The Lattice Landau gauge QCD 

We measured Kugo-Ojima confinement parameterfQ] of various configurations in Landau 
gauge[|], ||, [70. The Kugo-Ojima parameter is defined by the two point function of the covariant 
derivative of the ghost and the commutator of the antighost and gauge field 




(2.1) 

as c = — m(0). It is a scalar function at vanishing momentum in the continuum theory, but in the 
lattice simulation, we measured the magnitude of the right hand side of the eq.([0|) with \i = v 
polarizations of A v and D^, and then observed in the case of asymmetric lattices there is a strong 
positive correlation between the magnitude and the lattice size of the axis whose directions are 
perpendicular to the polarization. In Figure 1 the temperature and polarization dependence of the 
Kugo-Ojima parameter of MILC/y is presented. 

Data in Table 1 show that the parameter c of MILC2 / and MILC3 / are consistent with 1 while 
the data of quenched 56 4 configuration are about 0.8. Using the covariant derivative 

D^U X}fl )=S(U X}fl )d fl + [A Xtli ,-} 
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and the parameter e denned as 



e 

h = c and dimension d = 4, it follows that h = in the continuum limit means the horizon 

d 

condition^]. 
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Table 1: The Kugo-Ojima parameter of the quenched 56 4 lattice and that of the MILC2/ and MILC3/ . c x is 
the polarization along the spatial directions, c t is that along the time direction, c is the weighted average of 
c x and Ct, e/d is the trace divided by the dimension and h is the horizon function deviation. 

The running coupling in the MOM scheme in Landau gauge is given by the product of the 
gluon dressing function and the ghost dressing function squared: 

a s {q)=q 6 D G {q) 2 D A (q). 

We fit the scale by comparing with the perturbative QCD(pQCD) result at the high momentum 
region. 
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Figure 1: Kugo-Ojima parameter m(0) of 
MILC/ r configurations of T/T c = 1.02(blue di- 
amonds), T /T c = 1.23(red stars) and T/T c = 
1 .32(green triangles). 
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Figure 2: The running coupling a s (q)/n of 
MILC/ in Landau gauge The pQCD result of 
Nf = 3 (upper dash-dotted line) and Nf = 2 
(lower dashed line) and the extraction of JLab 
are also plotted(blue boxes). 



Recently an experimental extraction of an effective strong coupling constant was published 
from the Thomas Jefferson National Accelerator Facility(JLab) collaboration^]. Using the Bjorken 
sum rule at momentum, they predicted that the strong coupling constant approaches it. This pre- 
diction and data for Q > 0.4GeV shown in Figure 2 are consistent with our lattice results. The 
infrared suppression in the Landau gauge does not agree with JLab result. 
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3. The Lattice Coulomb gauge QCD 

We adopt the minimizing function F v [g] = \\A g \ \ 2 = J^-tr (a^ x /A 8 x •) , and solve dfAfat 



using the Newton method. We obtain e = — \ — d,A ( - from the eq. d,A,- + djDj(A)e = 0. Putting 

—oD 

g(x,t) = e e in U g i(x,t) = g(x,t)Ut(x,t)g^ {x + i,t) we set the ending condition of the gauge fixing 
as the maximum of the divergence of the gauge field over N 2 — 1 color and the volume is less than 
10~ 4 , Max x> a(divA) a (x) < 10 . This condition yields in most samples 

±£(<//vA^) 2 ~Kr 13 . 

a,x 

We measure the color-Coulomb potential by first defining 

vg ul (x,y) = ((a,x\^-\-d 2 )^- l \b, y )), 

where the outermost bracket denotes ensemple average, a,b are color indices, x = (x,t) and y = 
(y,t). We perform the Fourier transform 



xy 

= (K^l^)^ = (l^*{x)^(x)^ (3.1) 
where we define 

^£1^^)^ = 16,9) (3-2) 

h b *(y) = dpJT 1 \b,q) (3.3) 
We evaluate y b (y) as a solution of the differential equation for a plane wave source p(y) 

(-dD) V b (y)=p b (y) (3.4) 



and 



and define <j)^ b (y) = d^ b {y). Following the DSEQ, we express the color-Coulomb potential as 

Vcou, (q) = D G (q 2 )q 2 f(q)D G (q) (3.5) 

where Dc{q) is the ghost propagator and f(q) is the Coulomb form factor. 

The color-Coulomb form factors of MILC2/, MILC3/ and DWF m = 0.01 shown in Figure 3 
indicate that (Xcoul (g)depends on the kind of fermion. In the case of DWF, we excluded exceptional 
samples whose VqouI (9) at zero momentum is extremely large ( there are a few such samples in 
DWF m=0.01 among 50 samples). The results of MILCyy shown in Figure 4 indicate its tempera- 
ture dependence. 

The gluon propagator obtained as an average of data measured at each time slice does not 
depend much on the kind of fermion and temperature as shown in Figure 5 (MILC2/ and MILC3/) 
and in Figure 6(DWF). 
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Figure 3: The color-Coulomb form factor f(q) 
of DWF m = 0.01(green triangles), that of 
MILC2/(violet stars) and MILC3/(magenta dia- 
monds). 
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Figure 5: The Coulomb gauge transverse gluon 
propagator D^(q) of MILC2/ (violet stars) and 
M1LC3 /-(magenta diamonds). 




Figure 4: Same as Figure 3 but that of MILC/,, 
T/T c = 1.02(blue diamonds), T /T c = 1.23 (red 
stars)AC T/T c = 1.32(green triangles). 
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Figure 6: The Coulomb gauge transverse gluon 
propagator D&{q) of DWFm = 0.01 (green tri- 
angles), m = 0.02(magenta diamonds) and m = 
0.03 (orange stars). 



In the pQCD region the running coupling is defined from the color-Coulomb potential CCcoul (l) - 
1 lN \2N Nf 1 2 Vcoui (l) ■ In the infrared region, we measure the running coupling in Coulomb gauge as, 

a I (q)=q 5 D G (q) 2 D A (q) 



where Da (q) is the gluon propagator. 

The gauge field A Xjfl and the link variable U^(x) are related by, U^(x) = e "/< in [14], but 

JO 

Ufi(x) = e A X 'V- in our case, where A a = —= with Gell-Mann's SU(3) generator X. 

v2 

Conformance with the unit length transporter e~ A ^^ in the continuum theory yields the gluon 
dressing function DA{q) that approaches 2 in large q. Conformance with the transporter 
i.e. the normalization of A v ^ such that Da (9) approaches 1 yields the dj(q) consistent with the 
JLab result. We multiply 1/2 to Da (9) in the calculation of GCj(q). 

The running coupling aj (q) of MILC2 / and MILC3 / are shown in Figure 7 and that of MILC f t 
is shown in Figure 8. 

We observe that «/ (q) of MILC 2 /, MILC3 / 2 r 

and DWF (Figure 9) agree with the extraction of 
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Figure 7: The running coupling <Xi(q)/n of 
MILC3f and M1LC2/ in Coulomb gauge. The 
pQCD result of Nf = 3 (upper dash-dotted line) 
and Nf = 2 (lower dashed line) are also plotted. 
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Figure 8: The running coupling aj(q 2 )/n of 
MILC/,. r/r t = 1.02(blue diamonds), T/T c = 
1.23(red stars) and T /T c — 1.32(green triangles). 



JLab. We emphasize that the adjusted parameter 
is only the normalization of the gauge field which 
appear in the transporter as a coefficient of A. 

In case of zero temperature quenched configuration, the Kugo-Ojima color confinement pa- 
rameter c is about 0.8. In case of unquenched configurations it is consistent with 1 at zero temper- 
ature and decreases as temperature rises. 

The A 2 condensate is finite and ghost condesate parameters is small at zero temperature, but 
near T c they are consistent with 0. The gluon propagator in the Landau gauge and in the Coulomb 
gauge are infrared finite. Its large volume limit is yet to be investigated. 

The running coupling of the Coulomb gauge (Xi{q) may approach % at zero momentum, inde- 
pendent of the kind of fermion. 

The infrared suppression of the running coupling in Landau gauge may be due to the singular- 
ity of the color anti-symmetric ghost propagator that disturbs the color diagonal ghost propagator. 
Quark has the effect of magnifying the square norm of the color anti-symmetric ghost propagator 
and reduces its fluctuationflTl]]. We observed that the color diagonal ghost propagator is tempera- 
ture independent. Color anti-symmetric ghost propagator of a system containing dynamical quarks 
is temperature dependentjjT^]. 

Presence of dynamical quarks is crucial in defining infrared features of the QCD. Lattice data 
always suffer from finite size effect. Comparison with theoretical approaches as Dyson-Schwinger 
equation[]T0|] is important. Calculation of the effective coupling in other scheme and comparison 



with experimental data are also challenging issues [15|, |1 



The numerical simulation was performed at KEK using Hitachi-SRl 1000, at Yukawa institute 
of theoretical physics of Kyoto University using NEC-SX8 and at CMC of Osaka University using 
NEC-SX8. This work was supported by the KEK supercomputing project No.5(FY2006). We 
thank Alexandre Deur for sending the results of JLab, the MILC collaboration for supplying very 
useful data, and the BNL lattice QCD group for the information on the gauge configurations of the 
QCDOC collaboration. 
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